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Abstract 

A system of interacting, identical fermions described by standard Landau Fermi-liquid (FL) the- 
ory can experience a rearrangement of its Fermi surface if the correlations grow sufficiently strong, 
as occurs at a quantum critical point where the effective mass diverges. As yet, this phenomenon 
defies full understanding, but salient aspects of the non- Fermi-liquid (NFL) behavior observed be- 
yond the quantum critical point are still accessible within the general framework of the Landau 
quasiparticle picture. Self-consistent solutions of the coupled Landau equations for the quasiparti- 
cle momentum distribution n{p) and quasiparticle energy spectrum e{p) are shown to exist in two 
distinct classes, depending on coupling strength and on whether the quasiparticle interaction is 
regular or singular at zero momentum transfer. One class of solutions maintains the idempotency 
condition n?'{p) = n{p) of standard FL theory at zero temperature T while adding pockets to the 
Fermi surface. The other solutions are characterized by a swelling of the Fermi surface and a 
flattening of the spectrum e{p) over a range of momenta in which the quasiparticle occupancies lie 
between and 1 even at T = 0. The latter, non-idempotent solution is revealed by analysis of a 
Poincare mapping associated with the fundamental Landau equation connecting n{p) and e{p) and 
validated by solution of a variational condition that yields the symmetry-preserving ground state. 
Paradoxically, this extraordinary solution carries the burden of a large temperature-dependent ex- 
cess entropy down to very low temperatures, threatening violation of the Nernst Theorem. It is 
argued that certain low-temperature phase transitions, notably those involving Cooper-pair for- 
mation, offer effective mechanisms for shedding the entropy excess. Available measurements in 
heavy-fermion compounds provide concrete support for such a scenario. 

PACS numbers: TLlO.Hf, 71.27.+a, 71.10.Ay 



I. INTRODUCTION 

Landau Fermi liquid (FL) theoryl^^ is universally recognized as a cornerstone of modern 
low-temperature condensed matter physics. This theory predicts that the magnetic suscep- 
tibility x(T) becomes independent of temperature T as T — )■ 0. It also predicts that in this 
regime the entropy S(T) varies linearly with T, implying the same behavior for the specific 
heat C(T) = TdS{T)/dT and thermal expansion coefficient a = I'^dl/dT oc —dS{T)/dP. 
These predictions are in excellent agreement with available experimental data on conven- 
tional Fermi liquids, notably three-dimensional (3D) liquid ^He and the electron liquid 
present in ordinary metals. 

However, beginning the mid- 1990 's experimental studies of many-fermion systems have 
yielded abundant evidence of the failure of the FL picture upon entry into the regime of 
strong correlations among the particles. Such non-Fermi-liquid (NFL) behavior was first 
revealed and explored in experiments^"- on films of liquid ^He and found subsequently in 
electron systems in solids, especially heavy-fermion compounds.—"— One of the striking 
features of the observed NFL behavior is that as the temperature T drops to zero, both 
the spin susceptibility x(T) and the Sommerfeld ratio ■y(T) = C(T)/T = dS(T)/dT diverge 
when the density (or doping) reaches a critical value, while the Sommerfeld- Wilson ratio 
Rsw = x{T)/'y{T) changes much more slowly. 

The corresponding point in the Lifshitz phase diagram showing the sequence of different 
phases that replace one another at T = under variation of a suitable control parameter 
(coupling strength, density, doping) is generally called the quantum critical point (QCP). 
If the QCP were the end point, from finite T values, of a line of transition points Tc{B) 
associated with some control parameter B, then according to the theory of second-order 
phase transitions, divergence of 7(T) and x(T) at T — )■ would be specified by different 
critical indexes. For example, in the Landau theory of second-order phase transitions, the 
critical index q characterizing the divergence of the Sommerfeld ratio is zero, while that 
for the divergence of the magnetic susceptibility x is 4/3. In Wilson theory, the critical 
indexes have almost the same values as in Landau theory.— However, experiment deviates 
strongly from such theoretical predictions. For example, the critical-index value <;^ = 0.38 has 
been determined^ for the compound YbRh2Si2, which belongs to a family of heavy-fermion 
metals comprehensively studied by F. Steglich and his collaborators.—"— 



Further, the posited second-order phase transition, attributed as usual to critical spin 
fluctuations, possesses unusual properties: the magnetic moment specifying the ordered state 
of the metal YbRh2Si2 is extremely small, ~ 10^^ fiB, as if the order parameter were hidden.— 
Furthermore, in dense ^He films, where the emergent NFL behavior was documented for the 
first time, experiment has not identified any related second-order phase transition. Thus the 
theory of second-order phase transitions is hardly relevant to explanation of the observed 
NFL behavior. 

Consider, on the other hand, that within FL theory the quantities x ^ind 7 are both 
proportional to the density of states N{0). In turn, A^(0) is a linear function of the effective 
mass M*, which specifies, via the dispersion de{p — )■ PF)/dp = Pf/M*, the spectrum e{p) 
of single-particle excitations near the Fermi surface. Thus, the NFL behavior emerging at 
the QCP can be explained within the broader framework of FL theory by attributing the 
common divergence of x and 7 to divergence of the effective mass M* at the QCP. 

The NFL behavior of the entropy S, which is the subject of the present article, is exhibited 
primarily in the divergence, at the QCP, of the thermal expansion coefficient a{T — )■ 0). The 
most challenging data have been obtained for a group of heavy-fermion metals.— For CeCoIus 
in particular, a{T — )• 0) is found to be almost independent of temperature and exceeds typical 
values for ordinary metals by a huge factor 10^-10^, implying a corresponding enhancement 
of the entropy itself (see Fig. [6]). 

Such anomalous behavior of the entropy S{T -^ 0) is documented not only in the electron 
liquid within solids but also in liquid ^He films.- The empirical evidence implies the existence 
of a group of Fermi systems that exhibit paradoxical behavior of the entropy at the lowest 
temperatures currently accessible to measurements, behavior in contradiction to the famous 
Nernst theorem requiring S{T = 0) = 0. Our objective is to elucidate this behavior. 



which was already predicted two decades ago in Ref. |21| and was recently rediscovered in 
a very different theoretical context, ^^ namely the finite-charge-density sector of conformal 
field theory (CFT) based on the AdS/CFT gravity /gauge duality. Briefly stated, the origin 
of the NFL behavior of the entropy lies in the flattening of the single-particle spectrum that 
occurs near the Fermi surface when the strength of correlations attains a critical level that 
triggers this phenomenon, which may be envisioned as a swelling of the Fermi surface. 



II. REPRISE OF THE STANDARD FL QUASIPARTICLE PICTURE 

To set the stage for discussion of the issue of entropy excess in strongly correlated Fermi 
systems, let us recall the elements of the standard Fermi-liquid quasiparticle picture estab- 
lished by Landau.-i^ At the heart of this picture is the postulate that there exists a one-to-one 
correspondence between the totality of real, decaying single-particle excitations of the actual 
Fermi system and a set of immortal interacting quasiparticles, whose number is equal to the 
given number A^ of real particles. This condition is written as 

f N 

Tt J n{p)dv= — = p, (1) 

where n{p) is the quasiparticle momentum distribution, p is the density, Tr implies summa- 
tion over spin and isospin variables, and dv = (ip/(27r)^ is a volume element in a momentum 
space of dimension D. In FL theory, all thermodynamic quantities, such as the ground state 
energy E, the entropy S etc., are treated as functionals of the quasiparticle momentum 
distribution ra(p), the entropy S being given by the standard combinatorial expression 

S{n) = -Tr f (n(p) lnr2(p) + (1 - n(p)) ln(l - n(p))) dv. (2) 

In the case of homogeneous matter addressed in our treatment, a conventional variational 
procedure based on Eq. ([2]) and involving restrictions that impose conservation of the particle 
number and energy, leads to the connection^ 

n(p)= [l + e^(P)/^]"\ (3) 

between the momentum distribution and the quasiparticle energy 

measured from the chemical potential p. Assuming the quasiparticle spectrum to be T- 
independent, one sees immediately that at T — )■ the quasiparticle momentum distribution 
n{p,T) coincides with the Fermi step npip) = 6{pf —p)- 

However, as emphasized by Landau,- the quasiparticle spectrum e(p) depends on T, a 
dependence that happens to be crucial for resolution of issues addressed in this article. The 
structure of the spectrum e{p) can be derived from the Landau equation^"- 

5 



This relation expresses the deviation of the quasiparticle group velocity v{p) from the bare 
particle velocity p/M (with M the bare particle mass), in terms of the Landau interaction 
function /(p, Pi). The latter function is commonly treated as a phenomenological input 
and specified by a set of phenomenological parameters associated with the harmonics of its 
Legendre polynomial expansion. The three relations ([T]), (^, and (^ then determine the 
momentum distribution, quasiparticle spectrum, and chemical potential self-consistently. 

Half a century ago when FL theory was created, solutions of such a complicated integro- 
differential equation as ([5]) were unavailable. Fortunately, in ordinary Fermi liquids such as 
bulk liquid ^He or the conduction electron system in familiar metals, correlations between 
particles are moderate, allowing thermodynamic properties of these systems to be adequately 
described without the aid of a computer. Indeed, in these systems, the Fermi velocity 
vf = v{pf) = Pf/M* maintains a positive value, implying that the effective mass M* is the 
single parameter that the characterizes solution 

e{p)=pF{p-pF)/M* (6) 

of Eq. ([5]) in the domain |e(p)| < T relevant to the thermodynamics of Fermi systems. 

According to Eq.([5]), the quasiparticle effective mass of a 3D Fermi system is expressed 

as follows^i^ 

M _ I ,_ PfM 

in terms of the first harmonic /i of the interaction function /(p, pi), where the dimensionless 
parameter F^ = JiPfM/tt'^ has been introduced. Since the spectrum ([6]) differs from that 
of an ideal Fermi gas merely by the numerical factor M*/M, the low-temperature thermo- 
dynamic properties of conventional Fermi liquids — for which M* remains finite — behave in 
just the same way as those of the ideal-gas systems apart from trivial scaling. In particular, 
the spin susceptibility x{T) and Sommerfeld ratio 7(T) = C(T)/T differ from their results 
for the corresponding ideal-gas system simply by the factor M*/M. Furthermore, upon 
inserting Eq. ([6]) into Eq. ([3]) we find that the T = quasiparticle momentum distribution 
n{p) coincides with the momentum distribution npip) = d{pF — p) of the ideal Fermi gas. 
Referring to Eq. ([2]), it is seen that the entropy of the interacting system satisfies S{T = 0) 
in obedience to the Nernst theorem. 



III. DIVERGENCE OF THE EFFECTIVE MASS AT THE QCP 

This time-honored FL quasiparticle picture, in which a real system of interacting fermions 
is treated as a "gas of interacting quasiparticles,"— worked flawlessly for over three decades, 
but as pointed out in the introduction, more recent experimental studies, especially of novel 
materials, have brought indisputable evidence for its failure when the correlations grow too 
strong. Within FL theory, the strength of correlations is measured by the magnitudes of lead- 
ing harmonics of the dimensionless interaction function F = fN{0) where N{0) = ppM* /ti"^ 
is the density of states. The Fermi system is considered to be strongly correlated when one 
(or both) of the leading Landau parameters F^ or Fi reaches a value close to 1. Since the ef- 
fective mass diverges at the QCP, both these parameters diverge as well, so the QCP regime 
does indeed belong to a region of the Lifshitz phase diagram occupied by strongly correlated 
Fermi systems. The proportionality to A^(0) of both the spin susceptibility x{T) and the 
Sommerfeld ratio 7(T) implies that they become divergent as the temperature T goes to 
zero and hence exhibit NFL behavior. At the same time, the QCP Sommerfeld- Wilson ratio 
Rsw = x(T)/'y(T) remains finite. Additionally, it should be noted that the divergence of 
the specific heat at the QCP means that the entropy S{T) must also begin to participate in 
NFL behavior. 

The emergence of the QCP is a conspicuous signature of the essential qualitative differ- 
ence between Fermi systems that are strongly correlated and those with weak or moderate 
correlations. Three-dimensional (3D) liquid '^He, for which M* remains finite at any den- 
sity, belongs to the class of systems with moderate correlations. On the other hand, its 2D 
counterpart has been assigned to the class of strongly correlated Fermi systems based on 
evidence for the existence of a QCP from experimental studies of dense ^He films, as cited 
above. 

The quest for a fundamental understanding of QCP phenomena persists as one of the most 
important and most controversial activities in low-temperature condensed matter physics. 
Here we refrain from repeated disputation^iiM with respect to numerous papers (see, for 
example, Ref. l25|) based on the idea that the FL quasiparticle picture (in its wider scope) 
ceases to be applicable due to vanishing of the quasiparticle weight z in the single-particle 
state, inherent at points of second-order phase transitions. Aside from such theoretical 
considerations, this proposal is unable to provide a consistent description of experimental 



data in the QCP region. Rather, we pursue a different approach developed within the 



broader framework of the original Landau quasiparticle picture, dating back to Refs. 
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explored in detail in Refs. |29|-|36|. sketched above and encapsulated in Eqs. dl]), (E]), and 
i^. The emergent theory behaves properly on both the sides of the QCP,— i^i^ facilitating 
a consistent account of the properties of strongly correlated Fermi systems, including the 
entropy paradox. 

IV. POINCARE MAPPING IN THE THEORY OF STRONGLY CORRELATED 
FERMI SYSTEMS 

The development of a quantitative theory of strongly correlated Fermi systems within the 
quasiparticle framework rests upon solution of the fundamental equation (|5]) of FL theory, 
applied in the general context where n{k, T = 0) need not be the Fermi step np{p) = 6{pf — 
p). Iteration procedures provide a common strategy for finding solutions of such nonlinear 
integro-differential equations, and we will see that they also pave the way to a resolution of 
the entropy paradox that has arisen. Indeed, the iterative operation introduced in nonlinear 
dynamics known as Poincare mapping, and familiar in the theory of turbulence,—*^ will 
prove instrumental to elucidation of the striking features inherent in solutions of Eq. ([5]) 
beyond the QCP. 

To be more specific, the discrete iterative map corresponding to Eq. (^ reads 

the iterate jj^^^^^ for the chemical potential being determined from the normalization con- 
dition ([1]). The index j = 0, 1,2, . . . counts the iterations (zeroth, first, second, . . .). The 
iterate n^^'^^^p) of the momentum distribution n{p) is generated by inserting the corre- 
sponding spectral iterate e^^^^\p) into the right side of Eq. ([3]), which at T = reduces to 
a Heaviside function n{p) = 6{—e{p)). 

In canonical Fermi liquids, for which the sign of the Fermi velocity vp determined by 
Eq. ([7]) is positive and the function e'^'^\p) has the single zero aX p = pp, the first iterate 
n^^\p) and all higher iterates for the distribution n{p) coincide with np(p) = 6{pp — p), 
this being a fixed point of the transformation. However, beyond the QCP the sign of the 
group velocity vp becomes negative, and the first iterate e^-^\p) = e(p; np) for the spectrum 
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already has two or three zeroes, implying the emergence of new pockets of the Fermi surface, 
or equivalently, new kinks in the momentum distribution n^^\p). 

What happens in higher iterations depends on the structure of the interaction function 
/. Usually this function is smooth, i.e. regular, in momentum space, although the bare 



interaction may have long-range character in coordinate space. In t 
occurring beyond the QCP were first uncovered theoretically in Ref. 
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lis case, transitions 
and later addressed 



in Ref. 



39|. An important property of all these solutions is that they have a multi-connected 



Fermi surface satisfying the relation 

n'^(p) =n{p), (9) 

which is satisfied trivially by n{p) = 6{pf — p) in standard FL theory, continues to hold, 
implying that S{T = 0) = in agreement with the Nernst theorem. Thus, the QCP heralds 
a sequence of topological rearrangements of the Landau state in the strong correlation regime 
through an avalanche of new pockets of the Fermi surface. 

However, this scenario based on a proliferation of new pockets is not universal. There 
exists a different scenario for topological rearrangement of the Landau state, ^^ which was 
suggested even earlier: at some finite critical value of the coupling constant, the number 
of bubbles or "Fermi pockets" becomes infinite. The next subsection will demonstrate how 
this more radical scenario emerges within a suitable iterative procedure. (Additional details 
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40l .) It is this scenario that is directly relevant to explanation of 



may be found in Refs. 
the entropy paradox. 



A. 2-cycles in Poincare mapping for systems with a singular interaction function 

/ 

A remarkable feature of equations considered in modeling turbulence is the doubling of 
periods of motion in route to dynamical chaos.— "^^ The nonlinear system corresponding to 
Eq. ([8]) can be analyzed within this context by associating an iteration step with a step 
in time. One's first instinct is to object that such a phenomenon should not occur when 
Poincare mapping is implemented, and assert that iterations must converge, since chaos in 
the classical sense cannot play a role in the ground states of Fermi liquids at T = 0, which are 
assumed to possess a unique, nondegenerate structure. And indeed, this general assertion 



seems to be validated by comments made in the preceding subsection. However, this is no 
longer the case if the effective interaction / turns out to have long range in coordinate space, 
so as to become singular at the origin in momentum space. 

For example, the singularity that exists in the transverse bare interaction 

rOf n\ Pi ■P2- (pi-q)(p2-q)/g^ f^^. 

r^(pi,P2,q,w = o) = -^ (10) 

between quarks in dense quark-gluon plasma is not regularized by polarization effects, in 
contrast to the well-known regularization (screening) of the longitudinal Coulomb interaction 
due to gauge invariance. In this example, Eq. ([5]) takes the forn>^ 

^<P^ ^ + g f In ^P^^dp,, (11) 



dp dp J p — Pi dpi 

where e^ ~ cp is the bare single-particle spectrum of light quarks, with c the velocity of 
light. 

The first iterate for the spectrum, evaluated from Eq. (TTTj) with the distribution np{p), 
has an infinite negative derivative de^"* /dp at the Fermi surface. The corresponding first 
iterate of the momentum distribution is n^^\p) = 9{x + xq) — 9{x) + 9{x — xq), where 
X = p/pf — 1- The next iteration step yields n^'^\p) = nF{p), so the standard FL structure 
of the momentum distribution is recovered. The nonlinear system enters a 2-cycle that is 
repeated indefinitely. 

The first two iterations of the mapping process are depicted in Fig. [1] The top-left panel 
of this figure shows the bare spectrum e^^'{p) = e°(p). The first iterate of the spectrum, 
e^^\p), appearing in the middle-left panel, is evaluated by folding the kernel \n{2pF/(p — pi)) 
with the Fermi-step nF{p), shown in the top-right panel. The spectrum e^^\p) possesses 
three zeroes: pi < pp, P2 = Pf, and p^ > pp, implying that the first iterate n^^\p), drawn in 
middle-right panel, describes a Fermi surface with three sheets. This distribution differs from 
the ordinary Fermi step only in the momentum interval —Xq < x < Xq. The next iterate, 
e^'^\p) (bottom-left panel), again has a single zero pf, and the corresponding momentum 
distribution n^'^\p) (bottom-right panel) coincides identically with npip)- 

This example is not unique in exhibiting 2-cycle terminal behavior. A 2-cycle Poincare 
mapping also arises in treating the well-known Nozieres model,— for which the interaction 
function / has the limited singular form /(q) = {27T)^g6{q) with g > 0. In this model, the 
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iterative maps (depicted in Fig. [2]) are generated from the equation 

e^-'+i) (p) + /i(J+i) = pV2M + gn^^^ (p) , (12) 

along with the normahzation condition ([T]) for n^^^^\p) = 9(—e^^^^\p)). Here, the odd 
iterates n^'^^^^\p) of the momentum distribution deviate from the np{p) in the interval 
—g/4:e% < X < g/4:e'p, but in even iterations, the Fermi step reappears intact. 

Numerical analysis demonstrates that similar 2-cycles arise when Poincare mapping based 
on Eq. ([8]) is implemented for other systems possessing an interaction function singular at 
fc — )■ 0. In all these cases, the emergence of a 2-cycle provides an unambiguous signal of 
the instability of the standard Landau state and all its discrete modifications involving any 
finite number of pockets of the Fermi surface. 

B. Modified Poincare mapping and ne^v insight from chaos theory 

Apparently, the occurrence of persistent 2-cycles in the iterative maps of Eq. ([H]) prevents 
us from finding self-consistent solutions of Eq. ([5]) beyond the QCP for a specific class of 
Fermi systems possessing singular effective interactions. It can be argued, however, that this 
failure is a consequence of the inadequacy of the iterative procedure employed, which works 
perfectly on the FL side of the QCP. Indeed, a refined procedure that mixes iterations does 
allow one to avoid the 2-cycle terminal behavior. Nevertheless, the improved procedure still 
fails to yield a solution: once again the iterations do not converge,— although the pattern 
of their evolution becomes more complicated and — as will be seen — both intriguing and 
suggestive. 

By way of illustration, let us consider a refined Poincare mapping for the Nozieres model, 
with the same coupling parameter g = 0.26^, as before. Choosing a mixing parameter (, the 
equation 

e^'+^\p) + /i(^'+^) = (1 - C) [e^'Hp) + /^^'^] 

+C[pV2M + gn^^\p)] (13) 

is used to generate the iterative maps shown in Fig. [3l Recovery of the ordinary Fermi 
distribution npip), which is an inherent feature of the standard iteration procedure at the 
each even iteration (see Fig. [2]), no longer occurs; indeed, the sequence of iterations fails to 
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converge, even to a limit cycle. The number of sheets remains three at the second iteration. 
At the third, however, seven sheets of the Fermi surface emerge, and the number of sheets 
continues to increase in successive iterations. The same occurs in dealing with dense quark- 
gluon plasma. It is the nonconvergence of iterative procedures that compelled the authors 
to stop searching for solutions of Eq. ( ITT]) . 



of Ref . 
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Nevertheless, this obstacle can be overcome. A key to arriving at an acceptable, self- 
consistent solution lies in the opportunity to treat the number j = 1, 2, 3, . . . of the iteration 
as a discrete time step.— In this sense, the sequence of panels in the left column of Fig. [3] 
shows the "temporal" evolution of the quasiparticle spectrum. At any time-step t, the 
single-particle energy e{p, t) drops off steadily as t goes to infinity, while the sign of e(p) 
switches unpredictably in a finite region Q of momentum space adjacent to the Fermi surface. 
These erratic changes of sign, depending in practice on fluctuations of voltage, temperature, 
humidity, etc., that affect the computer executing the iteration code, induce unpredictable 
indeterministic jumps of the occupation numbers n{p,t) between the two values and 1. 
On the other hand, the lack of convergence of the iteration process may be attributed to the 
presence of a kind of "quantum chaos" associated with the basic equation of FL theory for 
a certain class of interacting many-fermion systems. A deeper understanding of the chaotic 
evolution of the corresponding systems in the fictitious time t may well yield unique insights 
into the role of chaos theory in the description of quantum phenomena. 

Significantly, as t goes to infinity, the region of momentum space in which iterations of 
Eq. (jS]) fail to converge tends to a definite limit Q. Entropy is an expected byproduct of 
chaos. Accordingly, we tentatively adopt fi as a special entropy that is associated with a 
Fermi system under conditions in which iteration of Eq. (|5]) does not converge to a solution. 
This naive Ansatz for the special entropy can be refined as follows. A "time-averaged" 
single-particle energy l{p) is defined by the standard formula drawn from statistical physics, 

T J 

I{p) = hm i / e{p, t)dt ^ hm - ^ e^^\p), (14) 

■?=o 

and a corresponding time average 

1 /• 1 ^ . 

n{p) = lim - / n{p,t)dt = lim -J2^^^\p) 0-^) 

J=o 

is constructed for the momentum distribution n{p). Invoking Eq. (IT2|) . the relation between 
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the two time averages is given by 

e{p) = p^/2M - /i + gn{p) . (16) 

The mixing parameter ( appearing in Eq. (TT3|) cancels out in deriving this formula. 

Wherever the iterations converge, the average n{p) only takes values or 1, and then 
e{p) is a parabolic function of p coinciding with the true single-particle energy. On the 
other hand, as seen from Eq. f lT^ and verified by results shown in the left panels of Fig. 131 
the function e{p) vanishes identically in the domain Q where iterations of Eq. ([5]) do not 
converge, so that 

n{p) = n,{p) = ^~P'^'^^ , pen. (17) 

The symbol ra* [p) has been introduced to denote a smoothed momentum distribution, deter- 
mined by averaging iterates e{p, t) for the single-particle energy according to the prescription 

In summary, by executing (i) a modified Poincare mapping procedure based on Eq. (II 3p 
along with (ii) "time- averaging" of iterates for e{p) and n{p) in the manner of Eqs. (TT4|) and 
f lTSj) . a self-consistent solution of Eqs. (II]), (jS]), and ([5]) that obeys the Pauli principle can in 
fact be found in the domain Q where the sequence of iterates fails to converge. The solution 
so obtained is independent of the parameters specifying the refined iteration procedure. The 
boundaries Pi and pf of the momentum interval pi < pp < pf defining the domain of non- 
convergence are determined by the conditions n^,{p.i) = 1 and n^{pf) = 0. The quasiparticle 
momentum distribution n{p) corresponding to the new solution, hereafter written as n^{p), 
is given by 1 and 0, respectively, aX p <pi and p > pj, and by Eq. (fTTj) in between. 

This program is implemented similarly in investigations of other Fermi systems for which 
iteration of Eq. (^ does not converge to a solution. The central quantity of the iterative pro- 
cedure, namely the averaged single-particle energy e{p), is constructed by means of Eq. ( TT4l) . 
Since this quantity vanishes identically in the domain Q where iterative solution founders, 
we immediately obtain the closed equation 

e{p,n,) = 0, pen (18) 

for the smoothed, NFL component of the momentum distribution n^{p). Recalling the 
definition (j4]) of the quasiparticle energy e{p), Eq. flTSj) can be recast in the variational 
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foriii^i 

which is free of any trace of "time-averaging." 

We can quantitatively characterize the excess entropy 5* of the presumed ground-state 
solution at T = by inserting the momentum distribution n*(p) determined from Eq. ( 1T9|l 
into the combinatorial expression (j2]). By the definition of n^,(p), the integrand in Eq. ([2]) 
vanishes outside the domain Q, so S^: must be essentially proportional to the volume of this 
domain (also called Q). Setting n^{p) inside Q equal to the typical value 1/2, one does 
indeed arrive at 5* = 20 In 2. The existence of such a nonzero entropy excess means that 
the rearranged T = ground state of the system corresponding to this solution of Eq. ([5]) is 
statistically degenerate. In the next section, it will be demonstrated that such a degeneracy 
inevitably emerges at some critical coupling strength in some region of the Lifshitz phase 
diagram for any type of interaction function /, whether or not the effective forces between 
quasiparticles are of long range. These findings may open a new thread in discussion of the 
elusive concept of "quantum chaos," which has often reduced merely to analysis of putative 
chaotic signatures in the spectrum of solutions of quantum one-body problems. 

V. FLATTENING OF THE SINGLE-PARTICLE SPECTRUM AS THE SOURCE 
OF EXCESS ENTROPY 

As seen above, the NFL results (e.g., Eqs. ( [T6|) and ( [T7|) ) obtained from application of the 
refined iteration procedure for solving the basic PL equation ([5]) can be interpreted as coming 
from solution of the variational equation for minimizing the Landau energy functional E{n), 
with the chemical potential fi serving as a Langrange multiplier to conserve particle number. 

A. Breakdown of the necessary stability condition for the Landau state and the 
emergence of flat bands 

We may now describe, within in a unified theoretical framework, how two alternative 
and qualitatively distinct topological transformations of the Fermi surface can take place 
when conventional FL theory, characterized by the zero-temperature Fermi-gas quasiparticle 
momentum distribution np{p) = 9{pf ~p), gives way to a more general theory of correlated 
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Fermi systems that maintains the quasiparticle picture. Which alternative is reahzed de- 
pends on the behavior of the interaction function /. Before reaching the QCP, the standard 
ground-state FL distribution np{p) apphes because the necessary stabihty condition (NSC) 
for the Landau state is satisfied. This condition reads simply^"^ 

SEo = I e{p- nF)Sn{p)dv > 0. (20) 

In words, the NSC is satisfied provided that for any admissible, particle-number conserving 
variation 6n{p) from the standard Landau quasiparticle momentum distribution np{p), the 
change of the ground-state energy Eq turns out to be positive. This is the case in naturally 
occurring Fermi liquids, such as 3D liquid '^He or low-density neutron matter, where one has 
e{p) = vf{p — Pf) with a positive value of the Fermi velocity vp = pp/M*. In this case, the 
NSC is always met and the Fermi surface is singly connected, since the sign of the variation 
6n coincides with sign of the difference p — Pf and hence with the sign of e(p). 

However, consideration of the full Lifshitz phase diagram anticipates situations in which 
at a threshold value of some input parameter, such as the density or magnetic field, there is 
a discrete change in the number of roots of fundamental equation^ 

e(p) = 0. (21) 

This change occurs, for example, at a critical density poo where the Fermi velocity vf{Poo) 
vanishes, or, equivalently, the effective mass M*(poo) diverges — i.e., at the QCP. Beyond this 
point, the number of roots of Eq. fl2T]) increases stepwise, implying that the Fermi surface 
becomes multi-connected. In other words, the minimum of E{n) leaves the Landau point 
UfIp) and moves to other boundary points of the manifold {n{p)} of admissible momentum 
distributions, at which n{p) continues to take only the two values and 1. In this situation 
FL theory still holds, at least when the temperature remains sufficiently low. This is the 
alternative that applies specifically to a regular interaction function /. 

A qualitatively different transformation appears when the coupling strength of a regular 
interaction, here denoted by A, is further increased to another critical value Ac, at which 
the minimum of E{n) at a given particle number eventually sinks into the interior of the 
manifold {n}. The true ground-state momentum distribution n{p) within the extended 
Landau quasiparticle picture then becomes a continuous function n^{p) of p determined 
by Eq. (IT9l) . the chemical potential p being identified as a Lagrange multiplier associated 
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with particle number-conservation. It will be convenient to refer to such distributions n*(p) 
as non-idempotent solutions of the variational problem flT9|) on the manifold {n}, meaning 
that they visit interior points < n{p) < 1 of the closed interval [0, 1], whereas boundary 
solutions n{p) have values restricted to or 1 and hence are idempotent in the sense that 
n^ = n. For a singular interaction, a nonidempotent solution appears already at the QCP. 
To illuminate the nature and conceptual status of these novel NFL solutions, we invoke 
a mathematical correspondence of the functional E{n) with the energy functional E{p) 
of statistical physics as taught from textbooks. If the interactions are weak, the latter 
functional attains its minimum value at a density p determined by the size of the vessel that 
contains it, and which it fills uniformly. In such cases, solutions of the variational problem 
evidently describe gases. On the other hand, if the interactions between the particles are 
sufficiently strong, there emerge nontrivial solutions of the variational condition 

6p{r) 

that describe liquids, whose density is practically independent of boundary conditions. 

The energy functional E[n\ of our quantum many-body problem must have two analogous 
types of solutions, with an essential difference: solutions n*(p) of the variational condition 
flT9|) must satisfy the condition < n{p) < 1 imposed by fermion statistics. This condition 
cannot be met in weakly correlated Fermi systems, but it can be satisfied in systems with 
sufficiently strong correlations. We illustrate these points by comparing the behavior of 
two models, one having a singular interaction functional / and the second having a regular 
interaction, smooth in momentum space.-^ Results based on the condition (TT9|) are presented 
in the upper panel (singular interaction) and lower panel (regular interaction) of Fig. |H 

For both models, the energy functional takes the form 

Ein) = J ^n{p) dv + -J /(pi - P2)n{pi)n{p2)dvi dv2, (23) 

the interaction function being given by 

/(Pi - P2) = I ^ I (24) 

IP1-P2I 

in the singular case (ffist model) and being represented by 

/(Pi - P2) = 7 Kt—^ (25) 

(Pl - P2)^ + K^ 
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in the regular case (second model), with k = 0.07 pp- 
For the first model, the variational condition (IT9|) reads 

P^ , X f ^(Pi) 



+ ^ / I : dvi = /i, (26) 



2M J |pi -p| 
as is consistent with equation ([5]). For the second model, this condition has the form 

"' +a/ , "(P;\ , dv, = ,. (27) 



2M 7 (Pi-p)2 + fi: 

The first of these equations is solved analytically by exploiting the aforementioned math- 
ematical correspondence between the functionals E{n) and E{p), where p{r) is the number 
density of a system of charged particles moving in an external potential field and interact- 
ing with each other via Coulomb forces. The result is simply n^{p) = 1/A for p < Pf and 
n=K (p) = otherwise, with the boundary momentum pf being determined by particle-number 
conservation, i.e., n*pj/37r^ = p. This qualifies as a legitimate non-idempotent solution of 
the quantum problem under consideration if it obeys the Pauli exclusion principle, which 
requires that the coupling constant A exceeds the critical value Ac = Aqcp = 1 coincident 
with the value of A at which the FL effective mass diverges. Evidently this solution does 
not apply to weakly interacting systems characterized by A < Ac, for which n{p) = np{p) is 
the proper solution. 

Turning to the second model, again there exists a non-idempotent solution of Eq. ( !27|) 
for every A value. However, as in the first model, this solution fails to meet the Pauli 
principle at small coupling constants, implying that in this case, the Fermi step riplp) must 
be the ground-state momentum distribution of the problem. However, in contrast to the 
first model, when A exceeds Xqcp (where the effective mass diverges and hence the Fermi 
step npip) can no longer be the true ground-state quasiparticle momentum distribution), 
the correct solution of the problem continues to reside on the boundary of the manifold {n}, 
implying a multi-sheet Fermi surface. This is so because the non-idempotent solution of the 
variational condition ( TT9l) still contradicts the Pauli restriction, and only when A attains a 
new critical value Acrit = 1.52 Aqcp where the Pauli principle becomes immaterial, does the 



n^, solution apply. (For more detail, see Refs. |29|J34J ). 
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B. Insights from introduction of an external magnetic field 

Imposition of a magnetic field H on systems without a spontaneous magnetic moment 
splits the quasiparticle system into two subsystems, whose momentum distributions n^{p) 
and n^{p) correspond to the ±1/2 projections onto the direction of the magnetic field. There 
is a consequent doubling of equilibrium condition ( !T9|) , implying 

6E{n+,n.) ^ 6E{n+,n_) ^ 
57i^{p) 6n-{p) 

Let us explore this extension in the framework of a modified Nozieres model^ in which 
both components / and g of the interaction function have the simplest form, /(pi — P2) oc 
g{pi — P2) oc (5(pi — P2), such that the set of equations (128|) becomes algebraic. At T = 0, 
one has 

2^^^ + 2M + ^-^+b) + ^-^-(p) = ^^ , 

-\i^bH +^ + ^nM + ^n^ip) = /i • (29) 

In writing these equations we assume sign of the spin-spin component to be repulsive. Oth- 
erwise the system acquires a spontaneous magnetic moment, and the corresponding analysis 
lies beyond the scope of this article. 

A salient feature of this problem is revealed by subtracting the lower equation of the set 
(!29|) from the upper one to obtain 

= fXBH + gin+{p)-n^ip)). (30) 

Since the right side of this equation is just the effective magnetic field acting on the quasipar- 
ticle in the medium, we see that the magnetic field is completely screened in the interior of 
the system, prompting an analogy with superconducting systems (explained as fermion-pair 
condensates). The solution 

fi-£/2M ubH 
n±{p) = 7 ^^~ ^^ 

of the set of equations 0291) exists in the momentum interval pf < p < pf where the Pauli 
restrictions < n±{p) < 1 are met. The expulsion of the external magnetic field remains 
in effect until its magnitude becomes so large that all of the occupation numbers (1311) reach 
limits enforced by the Pauli restrictions. With this stipulation, the phenomenon of complete 
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suppression of the magnetic field acting on a quasiparticle, while derived here within a special 
model, remains generally applicable — as demonstrated by numerical calculations based on 
the set (12^ . However, at finite temperatures complete suppression occurs only in fields that 
are sufficiently strong, i.e., such that ^bH > T; otherwise the screening effect disappears.— 

C. Fermion condensation and the Lifshitz phase diagram 

We now address more directly the properties of the NFL solutions of Eq. fll9p . Recalling 
the definition (^ of the quasiparticle energy, Eq. f lT^ is recast economically as 

e{p) = 0, pen. (32) 

This means that the spectrum e{p) of the strongly correlated Fermi system, which coincides 
with the solution of the problem obtained within the refined iterative procedure, must be 
completely flat in some domain Q of momentum space, implying that the Fermi actually 
swells to become a volume in momentum space. 

The set of states for which Eq. (13 2 p is satisfied has been called the fermion condensate 
(FC). The phase transition in which the Fermi surface swells from a line to a surface in 2D, or 
from a surface to a volume in 3D (discovered 20 years agcM-^'—) was originally called fermion 
condensation — and more recently, swelling the Fermi surface or emergence of flat bands. The 
rearrangements of the standard Landau state that occur at and following the QCP, namely 
(i) the generation of additional Fermi surfaces (or pockets) and (ii) fermion condensation 
(or formation of a fiat band) have been identified as symmetry-preserving topological phase 
transitions (TPT's).— i^ In finite Fermi systems, the phenomenon analogous to fermion 
condensation expresses itself as the merging of neighboring single-particle levels.— Interest 
in flat bands was ignited recently^""— in connection with a related opportunity to raise the 
critical temperature Tc for termination of superconductivity of the electron systems in solids. 

Another salient feature concomitant with flattening of the single-particle spectrum is 
partial occupation of single-particle states of given spin. As already seen, in the presence of 
strong correlations, the momentum distribution n{p) is no longer restricted to the values 
and 1, but may take any value within the interval [0, 1]. Aided by the variational condition 
( !T9|) . one can acertain what happens in this case when a quasiparticle with momentum p G f2 
is added to the system, again assumed to be homogeneous. In contrast to what occurs in 
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a canonical Fermi liquid, addition of just one quasiparticle now induces a rearrangement of 
the whole distribution function n^,{p) within the domain Q. It follows that the system can no 
longer be described as a gas of interacting quasiparticles, even though the original Landau 
quasiparticle concept still applies. 

In the intermediate region of the Lifshitz phase diagram, corresponding to Xqcp < A < 
Xpc in the second example of Section. V.A, the Fermi surface still remains multi-connected. 
A striking peculiarity of the segment of the Lifshitz phase diagram just beyond the QCP, 



revealed in Ref. |29j, is the fast breeding of new small pockets of the Fermi surface as A 
increases beyond Xqcp- The occurrence of such cascades is associated with the condition 
n{p) < 1 enforced by the Pauli principle. As we have seen earlier, solutions of the equilibrium 
equation flT^ always exist, but these solutions violate the restriction n{p) < 1 until the 
coupling constant A attains a critical value Xpc- It is this restriction that triggers a cascade 
of topological phase transitions in the segment of the Lifshitz phase diagram defined by 
Xqcp < X < Xpc- 

We are now ready to comment on the generic structure of the T = Lifshitz phase 
diagram of a homogeneous Fermi system. This diagram contains a FL domain, where 
correlations between particles are not so strong and the Fermi surface is singly connected. 
The quantum critical point separates this domain from a region where the Fermi surface 
acquires a new small pocket. In turn, this region adjoins a stretch of the diagram in which 
the Fermi surface contains two new pockets, and so on. In all these states, the FL relation 
11? {p) = n{p) of idempotency is preserved. 

The discrete generation of new pockets (or sheets) of the Fermi surface continues up 
to A = Xpc- At that point there emerges a new ground state that solves the variational 
condition flT9l) with a non-idempotent momentum distribution n*(p) satisfying the Pauli 
restriction < ^^(p) < 1 everywhere. This state wins the energetic contest with any state 
having a discretely multi-connected Fermi surface. It is this exceptional state that possesses 
a finite entropy excess 5**, in apparent contradiction to the Nernst theorem. 

So far we have restricted attention to the case of zero temperature. At finite T the 
degeneracy of the FC is lifted as an integral feature of the physics of strongly correlated 
Fermi systems,— and the spectrum e(p, T) acquires a finite slope with respect to T in the 
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domain Q. To determine this slope, Eq. ([3]) can be inverted to obtain 

.(p)=Tlni^. (33) 

Since a minute temperature elevation cannot affect the FC momentum distribution n^,{p), 
it can be inserted for n{p) on the right-hand side of this equation, yielding the dispersion 

of the FC spectrum at T — t- 0.^- The degeneracy of the FC spectrum is thereby removed, 
the FC dispersion now being proportional to T. 

A remarkable feature of the enlarged picture of the phase diagram, considered above 
at T > in the context of the FC solution, is that the behavior fl34p prevails not only at 
A ~ Ac, but already makes its appearance at A > Xqcp, i-e. just beyond the quantum critical 
point.— To make this point more quantitative, consider that in a rearranged ground state at 
T = having a multi-connected Fermi surface, the spectrum e{p,T = 0) shows drastically 
different behavior in different regions of the momentum variable p. It varies smoothly with 
p in regions away from the pockets of the Fermi surface, but oscillates very rapidly where 
the pockets reside, since the number of zeros of e(p, T = 0) must equal the number of 
pockets. The departure of |e(p, T = 0)| from zero in the latter domain, denoted below by 
Tm, emerges as a new energy scale of the problem, around which a "melting" of the structure 
associated with the discretely multi-connected Fermi surface takes place. Indeed, as seen 
from the basic Landau formula ([3]), the distribution n{p, T) remains almost unchanged from 
n{p,T = 0) as long as T < T^. However, as the temperature T increases, kinks in n{p,T) 
versus p are gradually smeared, and at T > T^, the function n{p,T) becomes continuous 
and almost T-independent in the momentum region where its abrupt variations occur at 
T = Tm- Accordingly, Eq. fl34|) comes into play and the dispersion of the single-particle 
spectrum e{p) becomes proportional to T.— In effect, at T > T^ the temperature evolution 
of e{p,T) and n{p,T) becomes universal, coinciding with their evolution in the state having 
aFC. 

Since the occupation numbers n^, are T-independent in systems having a FC, the entropy 
S{T —7- 0) determined by Eq. ([2]) has a nonzero value, as already foretold by the refined iter- 
ation/time average procedure for solving the basic equation ([5]). To illustrate this important 
property. Fig. ([5]) presents results of numerical calculations of the entropy S{T = 0) for the 
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model energy functionals ( 123|) and ( 125|) that contradict the Nernst theorem, which requires 
S{T — 7- 0) = 0. The contradiction is not evaded by going to finite T, even though the 
degeneracy of the single-particle spectrum is lifted. Elimination of the degeneracy does not 
lead to recovery of the basic FL relation n'^{p) = n{p), because the quasiparticle momentum 
distribution n{p) remains essentially the same as it was at T = 0. 

VI. RESOLUTION OF THE ENTROPY PARADOX 

The first possible route to removal of the entropy excess 5**, which is appropriate to the 
part of the Lifshitz phase diagram where the T = Fermi surface is multi-connected (i.e., 
^QCP < A < Xpc), has been already analyzed above, although not identified as such. The 
entropy S^, is released by means of a crossover at T ~ T^ from the FC state to a state with 
a multi-connected Fermi surface (or a Fermi surface decorated with pockets). It is likely 
that such topological phase transitions have already been observed experimentally^ii^ in 
the heavy-fermion metals CeCoIus and YbIr2Si2. 

Alas, this prescription for shedding the excess entropy is unavailable in the strong- 
coupling limit, when all the quasiparticles reside in the FC region. In this case, release 
of the entropy excess as T — )■ can occur in a different way. At this point it is important 
to recognize that in deriving the basic formulas dl]), (E]), and (^, we tacitly assumed that 
any phase transition involved in entropy release take places without breaking any symmetry 
inherent in the Landau ground state. With this in mind, resolution of the entropy paradox 
remains an open problem from the theoretical standpoint, though it is surely decided in 
nature from among a number of possible routes for draining the entropy excess. Different 
phase transitions exist as candidates even among those that preserve homogeneity of the 
system. Of these, it is instructive to focus attention on superconducting phase transitions. 
This route was already considered in this context in the first article^^ devoted to fermion 
condensation, because the preservation of homogeneity and isotropy in the superconducting 
ground state and the attendant violation of gauge invariance. By virtue of Cooper pairing, 
the spectrum E{p) of single-particle excitations of the superconducting ground state has a 
gap A. It is this gap that ensures a vanishing entropy at T = in agreement with the 
Nernst theorem, as seen in the BCS relation 

SiT ^ 0) oc e-'^^o)/^. (35) 
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Thus, one reliable way to eliminate S"* in systems with the well-developed FC is associated 
with Cooper pairing. The superfiuid/superconducting phase transition takes place provided 
the particle-particle interaction is attractive. In liquid ^He, this interaction is certain to be 
attractive in channels with large orbital angular momenta L.^^ At any rate, in 3D liquid ^He, 
pairing takes place in the triplet spin channel with L = 1. On the other hand, in electron 
systems of solids, pairing typically occurs in singlet channels. 

Within the framework of BCS theory, the critical temperature Tc of the superconducting 
phase transition in the singlet channel is determined from the BCS gap equation 

tanh(E(pi)/2T) 

mpi) 

where V denotes the effective interaction between quasiparticles in the corresponding pairing 
channel, while E{p) = [e^ip] + A|] . As usual, the momentum dependence of the gap 
function is neglected. At T — )■ Tc the gap A^ vanishes, and we are left with a linear 
integral equation for the single unknown parameter, the critical temperature T^. In the 
strong-coupling limit, where at T = all the quasiparticles belong to the FC, the critical 
temperature Tc is expressed in terms of the FC density with the aid of the identity 

tanh(e(p)/2Tj 1 - 2n,{p) 



Aip,T) ^lvip,p.)'-^^^^^§ffIlAip,,T)d^ (36) 



tip) 2Tc\n[{l-n,{p))/n4p)y 

Straightforward algebra leads to 



(37) 



In the strong-coupling regime, where the FC density coincides with the density p, the critical 
temperature Tc may then be estimated from this formula as 

Tc ~ Vp, (39) 

where V is a constant measuring an averaged pairing interaction. Since A(T = 0) ^ 1.76Tc 
in BCS theory, we infer that the critical temperature associated with release of the excess 
entropy is close to T^, which turns out to be of order of the Fermi energy, except for the case 
of the extremely small pairing constant. 

This situation applies in heavy-fermion metals, which possess electron systems having 
narrow 4/-bands, where NFL behavior — including the presence of the entropy excess — is 
well documented. However, in these systems the FC density is rather small and then, as 
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discussed earlier, the route to recovery of the Nernst theorem via the mechanism of singlet 
BCS pairing is not unique. Furthermore, this route can be blocked by imposition of a 
modest external magnetic field, sufficient to terminate superconductivity. As we have seen, 
imposition of an external magnetic field affects the FC only weakly, while the effect of the 
field on superconductivity, ultimately its destruction, is far more pronounced. Without 
attending to cumbersome formulas, we know that the critical magnetic field H^ sufficient 
for termination of superconductivity at zero temperature is proportional to the gap value 
Ai(0) oc Tc- Thus we can expect that imposition of modest magnetic fields H c:^ He on a 
superconducting system that hosts a FC in its normal state, i.e. at T > T^, causes re-entry 
of the system into such non-superconducting FC-bearing states at temperatures lower than 
Tc, thus moving the critical temperature for the onset of the Nernst regime toward zero. 

Is there any experimental evidence supporting the validity of this proposal, which pre- 
dicts quite unexpected behavior of strongly correlated Fermi systems at low temperatures? 
Analysis of experimental datafi^ on the thermal expansion coefficient a of the heavy-fermion 
metal CeCoIns reveals just such behavior in this material, which has the extremely low value 
Tc = 2.3 K for the critical temperature of its superconducting phase transition, compared 
with the characteristic energy scale of order 10^ K for the electron system of heavy-fermion 

are displayed in Fig. |6l We observe first that under de- 



metals. Data taken from Ref. 
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creasing temperature with no external magnetic field present, a does in fact undergo a jump 
at critical temperature Tc, affirming that one is dealing with a superconducting second-order 
phase transition, with a(T) subsequently dropping off to at T — )• 0, all as in the standard 
BCS regime. Second, imposition of a magnetic field H {B in the figure) to 4 T and then 
4.58 T shifts the jump of a(T) toward T = 0, in accordance with the BCS trajectory of the 
line Tc{H). 

However, interpretation of further details of the thermal evolution of the coefficient <y{T) 
presents a special challenge: at T > Tc = 2.3K the value of a{H = 0) is T— independent and 
greatly enhanced compared with values typical for ordinary metals. The behavior of a{H) 
under variation of an imposed magnetic field defies explanation within standard theoretical 
approaches. As H increases and Tc{H) correspondingly declines, the value of the jump in a 
at T = Tc grows rapidly, ultimately leading to recovery of the anomalous values of a that 
are found at T > Tc in the absence of the magnetic field (see Fig. |6]). We suggest that this 
perplexing NFL behavior can be explained in a natural way. Let us assume that the behaviors 
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in H and T of the thermal expansion coefficient and the entropy are ahke. Further, suppose 
the system carries an entropy excess due to the presence of a small portion of FC — which, as 
we have seen, is quite insensitive to the magnitude of H but highly sensitive to the presence 
of pairing correlations. It then follows that at H > He, the normal-state thermal expansion 
coefficient begins to decline only at temperatures lower than 1 K, furnishing evidence for the 
suppression of the entropy excess in harmony with the quest for satisfaction of the Nernst 
theorem. We suspect that this decline, which occurs without any jumps, can be explained as 
a crossover (or other transition) from the state with a FC to a state with a mult i- connected 
Fermi surface, as discussed above. 

Based on these experimental facts and theoretical considerations, we predict that analo- 
gous anomalies should be observed in measurements of the low-temperature thermal expan- 
sion coefficient of another heavy-fermion compound, namely the P-type of the heavy-fermion 
metal YbIr2Si2, whose entropy S{T)/N attains values around In 2 already at extremely low 
temperatures T ~ IK.— 

In systems with a well-developed FC, there exist other mechanisms that operate as effec- 
tively as singlet pairing correlations. Among them are triplet pairing correlations and the 
BCS-BEC crossover. However, as a preliminary analysis shows, the corresponding critical 
temperatures may be even lower than that associated with singlet pairing. 

VII. CONCLUSION 

We have explored the entropy paradox encountered in a many-fermion system when 
the correlations between particles become large enough to produce a rearrangement of the 
standard Landau state. The paradox arises from the fact that the system retains a large, 
temperature-independent entropy excess down to extremely low temperatures, behavior that 
is in apparent contradiction to the Nernst theorem. 

We have demonstrated that this unorthodox behavior can be explained within the general 
framework of the Landau quasiparticle theory. We find that it stems from a swelling of the 
Fermi surface that occurs when the correlation strength reaches a critical value. Equivalently, 
there is a flattening of the single-particle spectrum near the Fermi surface, with concomitant 
smearing of the ground-state quasiparticle momentum distribution. This phenomenon has 
been examined along two lines: 
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(i) Analysis of a Poincare mapping corresponding to the fundamental Landau equation, 
which connects the single-particle spectrum with the quasiparticle momentum distri- 
bution through the amphtude of the quasiparticle interaction, and 

(ii) Solution of the variational condition for the minimum of the ground state energy with 
a subsidiary restriction to ensure particle-number conservation. 

We have investigated routes for releasing the entropy excess as the temperature falls to zero, 
thereby avoiding violation of the Nernst theorem. Focusing on scenarios in which the system 
remains homogeneous, we have found that phase transitions associated with Cooper pairing 
are most effective in suppression of the entropy excess. We have shown that this conclusion is 
consistent with experimental data on the low-temperature behavior of the thermal expansion 
coefficient of the heavy- fermion metal CeCoIus. 
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